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Abstract 

We have performed normalization of Hamiltonian in the generalized photogravitational restricted 
three body problem with Poynting-Robertson drag. In this problem we have taken bigger primary as 
source of radiation and smaller primary as an oblate spheroid. Wittaker method is used to transform 
the second order part of the Hamiltonian into the normal form. 
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1 Introduction 

The restricted three body problem describes the motion of an infinitesimal mass moving under the gravita- 
tional effect of the two finite masses, called primaries, which move in circular orbits around their centre of 
mass on account of their mutual attraction and the infinitesimal mass not influencing the motion of the pri- 
maries. The classical restricted three body problem is generalized to include the force of radiation pressure, 
the Poynting-Robertson effect and oblateness effect. 

J. H. Poynting (1903) considered the effect of the absorption and subsequent re-emission of sunlight by 
small isolated particles in the solar system. His work was later modified by H. P. Robertson (1937) who used 
a precise relativistic treatments of the first order in the ratio of he velocity of the particle to that of light. 

The effect of radiation pressure and P-R. drag in the restricted three body problem has been studied by 
Colombo et al. (1966), Chernikov Yu. A. (1970) and Schuerman (1980) who discussed the position as well 
as the stability of the Lagrangian equilibrium points when radiation pressure, P-R drag force are included. 
Murray C. D. (1994) systematically discussed the dynamical effect of general drag in the planar circular 
restricted three body problem, Liou J. C. et al. (1995) examined the effect of radiation pressure, P-R drag 
and solar wind drag in the restricted three body problem. 

Moser's conditions (1962), Arnold's theorem (1961) and Liapunov's theorem (1956) played a significant 
role in deciding the nonlinear stability of an equilibrium point. Applying Arnold's theorem (1961), Leontovic 
(1962) examined the nonlinear stability of triangular points. Moser gave some modifications in Arnold's 
theorem. Then Deprit and Deprit (1967) investigated the nonlinear stability of triangular points by applying 
Moser's modified version of Arnold's theorem (1961). 

Bhatnagar and Italian (1983) studied the effect of perturbations on the nonlinear stability of triangular 
points. Maciejewski and Gozdziewski (1991) described the normalization algorithms of Hamiltonian near 
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an equilibrium point. Nicdzielska (1994) investigated the nonlinear stability of the libration points in the 
photogravitational restricted three body problem. Mishra P. and Ishwar B.(1995) studied second order 
normalization in the generalized restricted problem of three bodies, smaller primary being an oblate spheroid. 
Ishwar B.(1997) studied nonlinear stability in the generalized restricted three body problem. 

In this paper normalization of Hamiltonian is performed in the generalized photogravitaional restricted 
three body problem with Poynting-Robertson drag. Wittaker method is used to transform the second order 
part of the Hamiltonian into the normal form. 



2 Location of Triangular Equilibrium Points 

Equations of motions are 

where, V.-"*-™ (1) 
ox rf 

y + 2nx = U y , U y = — - — 3— (2) 

n 2 (x 2 +y 2 ) (1 -//)gi M M2 

Ul = n 1 1 r TTT 

2 n r 2 2r 2 

r\ = (x + nf + y 2 , r 2 = {x + n-l) 2 + y 2 , n 2 = l + ^A 2 , 
J y 1= (» + A*)[(» + M)» + w] + -_ Ar2= ^ + ^ + ^] + . + 

^ _ (i-M)(i-gi) _ "2 < I rni,m,2 be the masses of the primaries, A 2 = r \„l p be the oblateness 
coefficient, r e andr p be the equatorial and polar radii respectively r be the distance between primaries, 
q = be the mass reduction factor expressed in terms of the particle's radius a, density p and radiation 



pressure efficiency factor \ (in the C.G.S. system) i.e., q = 1 — ' : * ap — -. Assumption q = constant is 



5.6xl0~°x 
ap 

equivalent to neglecting fluctuation in the beam of solar radiation and the effect of solar radiation, the effect of 
the planet's shadow, obviously q < 1. Triangular equilibrium points are given by U x = 0, U y — 0, z = 0, y 7^ 0, 
then we have 



a;, = x <!l ^ - y - \ (4) 

y *- y Y 3M1 - M)% 3 ^ J ( ) 

where xo = 4j — /U, yo = ±6(l — x)^ 2 an< ^ ^ = 9i^ 3 ' as m preprint Kushvah & Ishwar (2006) 



3 Normalization of Hamiltonian 

We used Whittakar (1965) mathod for the transformation of £f 2 into normal form The Lagrangian function 
of the problem can be written as 

L=l(x 2 + y 2 ) + n(xy-xy) + ^(x 2 +y 2 )+ { l-^ + ^ + ^ (6) 
+ ^{ (X + ^ + ^ -narctan- V ^ 



2r\ " (x + fji). 
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and the Hamiltonian is H = —L + P x x + P y y, where P x ,P y are the momenta coordinates given by 



dL , W lf . „ dL . W 1 



For simplicity we suppose q\ = 1 — e, with |e| << 1 then coordinates of triangular equilibrium points can be 
written in the form 

7 e A 2 A 2 e (9 + 7) „, 476 TIT 



V3r 2e A 2 2A 2 e (1 + 7) 4 7 e -i 

V = "rH 1 - ~z ^ ~ ^nWi j=nWi > 

2 1 9 3 9 9^3 27V3 J 



(8) 



where 7 = 1 — 2jU. We shift the origin to L 4 . For that, we change x ^ x* + x and y — > 2/* + y. Let 
a = x# -\- [i)b = y* so that 

a _ l/ 2e „. , 2A 2 e (9 + 7)„ w 8 7 e 



(10) 



2 ( 3 3 3^ 27V3 

V3r 2e A 2 2,4 2 6 (1 + 7) 4 7 e 

& = — n - h -=-nWi j=nWi \ 

2 I 9 3 9 9V3 27V3 > 

Expanding L in power series of x and y, we get 

L = L Q + L 1 +L 2 +L 3 + --- (11) 
H = H + H 1 +H 2 +H 3 + --- = -L + P x x + P y y (12) 

where Lq, L\, L 2l L 3 . . . are 

3 2e 7 e 3 7 A 2 3A 2 e 

Ln = 1 r yA 2 

2 3 3 4 2 ' 

VSnWi 2 7 neVFi 23enWi & 

1 r=nWi = = narctan- (13) 

4 3V3 3^ 54V3 a 

. r V3 5A 2 7eA 2 AnW\ ljnWi -, 

Li = x-i = H = H V 

1 2 8V3 12V3 9 18 ' 

.,1 e A 2 eA 2 nWi 2nenWi, 

+ V ^2 ~ 3 ~~8~ + 12\/3~ 6vf + 3^3 ' 

, 1 7 9A 2 157A 2 35eA 2 297e.4 2 

~ ^ ~ 2 + 2 + IT + 8 12 12 - 

ZVZnWi _ 2£L n w - 5enWl - 7 7 enT ^i i 
8 3\/3" 1 12V3 4^3 ' 

r 15V3A 2 9V3 7 A 2 /- * 43e , 

- y{ ^ + g 2V3eA 2 - 2V37eyl 2 — + 7 nWi - — nWij (14) 

L 2 = I ?/) + n(xy - xy) + ^-(x 2 + y 2 ) - £x 2 - Fy 2 - Gxy (15) 
L 3 = -i {^Tj + 3x 2 yT 2 + 3xy 2 T 3 + y 3 T 4 + 6T 5 } (16) 
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where 



16 1 2 6V3 27^3 

r„ (199+ 17 7 ) Trr 2(226 + 99 7 )e TTr , 1 



(17) 



F = 



-1 
"16 



f-,^ „ , 717A 2 e (67+ 19 7 ) Tlr 2(413 - 397)e 
1 10 - 2e + 21^ 2 TT7 ^- - ^ — ^ f; nVFi + v ^ " nW x 
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G = 



, 293A 2 e (187+27 7 ) TI7 
37^ 2 e (13 + 7) 



6\/3 
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27^3 



27^ 



^{2 e + 6A 2 - 
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2V3 
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2V3 



27^ 
(186 -7)e 



}} 



(18) 



(19) 



Ti = 



To = 



T 3 = 



n = 



16 c , 979 , (143 + 9 7 ) „, (459 + 376 7 ) „, 



ii, 4e or, 1507 , (215 + 29 7 ) „, 2(1174+ 1697) „, 1 
-7 U4 + y + 25^2 - ^A 2 e - ^-^n Wl - W l£j 



3V3 
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,,16 A 2 367 , 115(1 + 7) _ _ _ (959- 136 7 ) _ _ _ 
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9 2 54 
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2 "3 £ 



23A 2 (37 + 7) (219 + 253 7 ) 
— 44A 2 e - -p^-nWi - - = — —nWie 



18V3 
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I 27 18V3 8lV3 



{ax + by) {3{ax + by) - (bx - ay) 2 } - 2{xx + yy){ax + by)(a 2 + fe 2 ) 



2(a 2 + b 2 f 

The second order part H 2 of the corresponding Hamiltonian takes the form 



(20) 
(21) 

(22) 

(23) 

(24) 



H 2 = Px + p y + n (y Px - xp y ) + Ex 2 + F 2 + Gxy (25) 

To investigate the stability of the motion, as in Wittaker(1965), we consider the following set of linear 
equations in the variables x,y: 



XP X 



-XR, 



i.e. AX 



dH 2 

dx 

dH 2 
dy 





Ax - ^ 

AX ~ dP x 



Xv = 2»a 



(26) 
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X = 



y 

Per 



and A = 



2E G A -n 

G 2F n A 

-A n 1 

-n -A 1 



(27) 



Clearly \A\ = implies that the characteristic equation corresponding to Hamiltonian H 2 is given by 
A 4 + 2{E + F + n 2 )A 2 + 4EF - G 2 + n 4 - 2n 2 (E + F) = (28) 

This is characteristic equation whose discriminant is 

D = A{E + F + n 2 ) 2 - 4{4EF - G 2 + n 4 - 2n 2 (E + F)} (29) 

Stability is assured only when D > 0. i.e 

H < Mc - 0.221895916277307669e+ 2.1038871010983331A 2 

+ 0.493433373141671349eA 2 + 0.704139054372097028nVFi + 0.401 154273957540929neVKi 

where fi Co = 0.0385208965045513718 

When D > the roots ±iuu\ and ±iw 2 (^1,^2 being the long/short -periodic frequencies) are related to 
each other as 



2 , 2 

UJ 1 + U! 2 



2 2 
^1^2 
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16V3 



16 



288V3 



288V3 



48V3 



(30) 



(31) 



(0 < lo 2 < —1= < uji < 1) 
v2 



From ( 30)and ( 31) it may be noted that u)j(j = 1, 2) satisfy 
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4e 107eA 2 
+ ~9 27~ 
16 32e 
+ ' ~27 + 243 + 
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+ -- 
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486V3 



486V3 
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81 
107nV^i 



729 243 
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243^ 27V3 
16 32e 208^2 
27 



27 

15892enVFi 
2187V3 729V3 
1880eA 2 2720nVKi 49552enVKi 80 7 enVyi\ 4 
2187V3 6561V3 2187V3 J j 



243 81 729 
Alternatively it can also be seen that if u — ujilo^, then equation ( 31) gives 



(32) 



2 4e 107e,4 2 

7 = 1 H 

' 9 27 



25nWx 3809enVFi /2e 1579eA 2 55 7 nT4^i 4961 7 enT^i 
27V3 486V3 7 \~3~ 324 9^3 486^ 



+ 



16 32e 208A 2 1880e^ 2 320nl^i 15856enVKi\ 2 
27 + 243 + ^1 729 + 243V3 ~ 2187^3 J 



(33) 
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Following the method for reducing H 2 to the normal form, as in Whittaker(1965),use the transformation 



X = JT where X 



x 

y 

Px 

Py J 



' J — [Jij]l<i<j<4, T 



Qi 
Q 2 
Pi 
Pi 



(34) 



P, = (2fe) 1/2 cos^, Q, = (^£)i/2 sin 4j (i = l,2) 



(35) 



The transformation changes the second order part of the Hamiltonian into the normal form 

H2 = OJ\Ii — U2I2 

The general solution of the corresponding equations of motion are 



(36) 



const., 



±ivi + const, (i = 1,2) 



(37) 



If the oscillations about L4 are exactly linear, the equation( 37) represent the integrals of motion and the 
corresponding orbits will bi given by 



= J13 y / 2uj 1 I 1 cos 4>i + J 14 y/2uj 2 h cos 4> 2 



y = 




— sin^i + J 2 2\j — sin 02 + J 2 aV^h LJ i cos^i + J 2 4 ^2/2^2 sin cj> 2 
wi V ^2 



(38) 
(39) 



where 
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with Z? 



4^ + 9, (j 



1,2) and fcf = 2^ 



lj ^2 



-2w| + 1. 



4 Conclusion 

Using Wittakcr(1965) method wc find that the second order part H 2 of the Hamiltonian is transformed into 
the normal form H 2 = u>\I\ — oj 2 I 2 . 
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